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Abstract 

Recently the partial wave cutoff method was developed as a new calculational scheme for a 
functional determinant of quantum field theory in radial backgrounds. For the contribution given 
by an infinite sum of large partial waves, we derive explicitly radial WKB series in the angular 
momentum cutoff for d = 2,3,4 and 5 {d is the spacetime dimension), which has uniform validity 
irrespectively of any specific values assumed for other parameters. Utilizing this series, precision 
evaluation of the renormalized functional determinant is possible with a relatively small number 
of low partial wave contributions determined separately. We illustrate the power of this scheme in 
numerically exact evaluation of the prefactor (expressed as a functional determinant) in the case 
of the false vacuum decay of 4D scalar field theory. 
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I. INTRODUCTION 



Functional determinants of (ordinary or partial) differential operators arise in many areas 
of physics: for instance, in connection with the one-loop effective action in quantum field- 
theoretic studies and in the semiclassical approximation to quantum mechanical tunneling 
amplitudes. However, explicit evaluation of these quantities, especially with partial differ- 
ential operators involving nontrivial background fields, is usually a very difficult problem. 
Explicit analytic results are known only in some simple cases, such as the one-loop effects in 
constant electromagnetic fields in QED [l, 2, 3] and in a covariantly constant field strength 
in non-Abelian gauge theories ji, 5|. Therefore various methods for approximate calculation 
were considered, the large mass expansion [6] and the derivative expansion [7,8] being good 
examples of them. But the validity of these approximate methods crucially depends on the 
range of various parameters entering the problem. 

Recently there has been a significant progress in this problem, at least when the differ- 
ential operators are separable. Especially, for background fields having radial symmetry, a 



le partial wave analysis has been developed in the form of the partial wave 



method using t 
cutoff method 

determinant \^ for an arbitrary value of quark mass. (The same quantity with massless 



10|]. This method was first used in the computation of QCD instanton 



quarks was calculated in a classic paper by 'tHooft [IfJ long time ago). It was then applied 
to the evaluation of the one loop effective action for more general classes of radial background 
fields [1^. The prefactor in the false vacuum decay rate, which requires an evaluation of the 
functional determinant also, is calculated by the same method 

Of crucial importance in the above-mentioned calculational scheme is to find a simple way 
to extract a finite renormalized quantity from the infinite sum of partial wave contributions. 
In {9] this was achieved by introducing a cutoff in the partial wave sum and then finding a 
uniform radial WKB expansion for the sum of partial waves beyond the cutoff value (which 
is combined with the conventional renormalization counterterms) . Combining the leading 
terms of this WKB expansion with the contribution from partial waves below the cutoff 
value, it is possible to secure a finite renormalized value in the limit of large cutoff value. 
[In ij], similar results were obtained using the zeta function technique]. More recently it 
is observed that the inclusion of higher order terms in the uniform radial WKB expansion 



greatly improves the rate of convergence for the infinite sum of partial wave contributions 
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[l5l |. The computational labour needed for the functional determinant calculation is thus 
much reduced. Efficiency of this scheme will become especially conspicuous for functional 
determinants of higher-dimensional differential operators, thus making it an effective tool 
also for the studies of higher-dimensional quantum field theories. 

In the present paper we will give a simplified derivation of the uniform WKB expan- 
sion and provide explicitly several leading terms of this expansion (needed for fast precision 
evaluation of functional determinants) in general contexts. It is our hope that these ex- 
plicit results find useful applications in various related problems. This paper is organized 
as follows. In Section [III the partial wave cutoff method is explained briefly and the de- 
sired form of the asymptotic WKB series is presented. In Section llllt after introducing 
the proper-time representation of the radial functional determinant, we derive the large / 
expansion of the proper-time Green function. In Section HVl the infinite sum of contributions 
from high angular momentum is explicitly evaluated using the radial WKB series and the 
Euler-Maclaurin summation method and then the uniform WKB expansions (as descend- 
ing series in the angular momentum cutoff L) are presented in various dimensions i.e., for 
d = 2,3,4,5. Our formulas for the renormalized functional determinants have definitely 
faster convergence property compared, say, to those given in [l^. In the next section we 
consider a direct application of these results, finding a more accurate value for the false 
vacuum decay rate in the context of 4D scalar field theory. In Section |Vl] we consider the 
functional determinant in gauge theories, where the radial potential has a linear dependence 
in angular quantum number /. But this does not change the general structure, and in this 
case also the appropriate coefficient functions in the uniform WKB expansion can be found 
in explicit forms. 

II. SETTING UP THE PROBLEM 

In order to set the problem precisely, let us start by considering a pair of partial differential 
operators 

M = -d^ + V{r), M^''"' = (2.1) 

where d"^ = d^d^ is the Laplace operator in d dimension and V{r) is a radial potential 
vanishing sufficiently fast at infinity. In the one- dimensional case (i.e., with -M = —-£2 + 



3 



V{r)) with the Dirichlet boundary condition on the interval [0,oo), we can determine the 
ratio of two functional determinants of the operators with mass m, using the Gel'fand and 



Yaglom's theorem 16|], as 

det{M + m^) ^ ^(oo) ^2 2) 

det(A<fr^^ + m2) ^/'^^^(oo)' ^ ' ' 

where the wave function ilj{r) satisfies the ordinary differential equation (ODE) {A4+m'^)ip = 

with initial value conditions at r = 0: 1^(0) = and ip'{0) = 1. The other function ip^'^'^'^lr) 

is the solution to the differential equation [Ai^^'^^ + iv?)ip^^'^'^ = with the same initial 

conditions. This method turns the problem of finding an infinite number of eigenvalues into 

that of finding the solutions to the ODE initial value problems. 

Now we consider the case of higher dimensions (i.e., d > 2). [In this paper we will provide 

explicit formulas for the cases with d = 2,3,4,5 but the extension to higher dimension is also 

straightforward]. When the potential is radial, i.e., V = V{r), we can use the partial wave 

analysis, taking advantage of the spherical symmetry. Formally, for the radially separable 

operators given in (12. ip . the logarithm of the determinant ratio can be written as a sum of 

the logarithm of radial (that is, one-dimensional) determinant ratios: 

\det[M^'<'^ + m^] J f^Q Vdet[A1f^^ + m2]y ^ ^ 

Here I denotes the angular momentum quantum number appropriate to each partial wave 
and 

(2l + d-2)(l + d-3)\ 

= — w^^. — p-^' 

is the degeneracy factor [IJ, [l8|. The associated radial differential operator A^/ is given by 



= -^r^rV Q-A^ ,2 +^(^)- (2-5) 



and Alf*^^ has the same form as M.i but without the potential term Vir). 

The individual radial determinant ratio in (12.31) is finite and it can be evaluated easily 
by using the above Gel'fand- Yaglom method; but the sum over I leads to a divergent result. 
This problem is related to renormalization and an elegant method to extract the finite or 
renormalized expression l^ren from (after a suitable regularization and renormalization) 
is presented in jo], [l^. We are concerned here with more practical problem, which should 
be addressed if one wants the full, including the finite part, expression of Oren- The rate of 
convergence of the Z-sum in (12. 3p is quite slow, and we require an efficient method to deal 



with this /-sum. To this end it is convenient to introduce a partial wave cutoff L [lOj and to 
spht the sum into two pieces, i.e., the low angular momentum part and the high angular 
momentum part f^H^ 

firen = ^^L + ^^H, (2.6) 

where 5VL denotes the 'conventional' renormalization counterterm. Separate treatment of 
and constitute the crux of the partial wave cutoff method. 

The part VLi^ (see (12.71) ) can be evaluated with the help of the Gel'fand-Yaglom method. 
Since the determinant ratio for given / behaves like ~ y for large / and the degeneracy factor 
gi{d) increases as Z'^"^, it should be clear that I^l behaves like L'^"^ for d > 3 and like InL 
for d = 2 m the large L limit. (This reveals the divergent structures in the formal expression 
in (12.31) ). As for the part which involves the sum of all partial wave contributions with 
/ > L + 1, we can evaluate it analytically in a uniform asymptotic series of the form 



j-OO ( °° \ 

dr\ Qiog + ^ Q-n^"" > (2-9) 

•'^ \ n=2-d ) 

where Q-nS may have an implicit L dependency of 0{LP) and Q\og behaves as O(lnL) in 
the large L limit. This uniform nature makes also the r integrals in (12. 9p well-defined. To 
find explicit forms of the Q's, we take the proper-time representation for the functional 
determinant of radial operators for each partial wave and then use the quantum mechanical 
radial- WKB expansion which becomes exact in the large / limit. For the desired large L 
expansion we then perform the sum over / = L + 1, ■ ■ ■ , oo with the help of the Euler- 
Maclaurin method. These are given in following sections. 

Note that, as L oo, unsuppressed terms in the expansion (12.91) may grow like 
InL, L, ■ ■ ■ , but they match precisely the large-L divergences coming from VL\^ ex- 

cept for the sign. Hence, combining the large-L-unsuppressed terms of fin with f^L and 
taking L — > oo limit, we get a finite renormalized quantity l^ren, i-e., 

d~2 



^i-en = lim 

L— >oo 



[ dri Qiog + QnL 

•^0 V n=0 



(2.10) 
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In [l^, Dunne and Kirsten identified this expression for d = 2, 3, 4 by using tlie zeta function 
technique. 

In principle, since fl2.1UI) yields a well-defined expression, one can use this expression 
to obtain the renormalized functional determinant. But we still have a practical problem 
determining Qj^. Since it is generally not possible to find a master formula for the determi- 
nant ratio valid for all /, we have to evaluate (numerically) those partial wave contributions 
corresponding to the angular momentum range < Z < L, to be able to determine JIl. 
Because of the slow rate of convergence, a very large number of these determinant terms 
should be thus considered to get a sufficiently good result for the sum. There is a rather 
simple way to secure a reliable large-L limit value in fl2.10l) with a relatively small number of 
partial wave contributions. Including the —suppressed terms of the expansion (12.91) inside 
the squared braces in (I2.10p would make the sum converge faster, thereby reducing the num- 
ber of partial-wave determinants to be evaluated explicitly. To appreciate this better, note 
that if we were able to calculate both of I^l and Qu exactly, their sum would be independent 
of the choice of the cutoff L. This implies that, leaving aside possible numerical inaccuracy 
in calculating ^l, the L-dependency in the sum for finite value of the cutoff L is really due 
to our ignoring of the -^-suppressed contributions in the asymptotic series (12. 9p . Therefore, 
it is possible to reduce this undesired L-dependency systematically by taking into account 
the (ignored) higher order terms in the -^-series. Now, instead of taking the strict L —>■ oo 
limit in (I2.10p . we can write the following formula for Q^en'- 

d-2 N 



a + I dr [q,,, + X: Q„L" + E Q-n-^j + O (^jj^j , (2.11) 



where refers to the order of truncation. In this formula the error is indicated by the last 
term and it is totally under control. It is apparent that, for a given value of the cutoff L, 
we get more accurate value of firen by taking into account more ^^-suppressed terms. Or, 
for a given accuracy, we can lower the value of L by including some -^-suppressed terms. 
We can thus use (12.111) as the basis of precision calculation for functional determinants. In 
this work, we take = 4 for concreteness and derive the expressions for Qd-2, ■ ■ ■ , Q-4 and 
Q\og in dimensions d = 2, 3, 4, 5 to facilitate the use of our formula (12. lip in various physical 
problems. 
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III. THE PROPER-TIME RADIAL GREEN FUNCTION AND ITS DERIVATIVE 
EXPANSION 



First we write the partial-wave determinant ratio in a more convenient form 
where 



(3.1) 



Ml = r('^-i)/2>l,r-(^-i)/2 (3.2) 



+ ' "\ ' ' +V{t) (3.3) 

and Alf®*^ is equal to J^i with V = Note that the operators J^i and Jvif^^ do not involve 
any first order derivative term. It is also convenient to introduce the effective potential 

V,(r) = ^ "^^i^ + + l-(r). (3.4) 

The Schwinger proper-time representation for the determinant ratio, for a partial wave 
/, is given as 

det(A^f«^ + m2) Jo s Jo I ' ' ^ i v , . ;/ v ; 

Here the proper-time radial Green function is defined by 

Ai{r,r';s) = {r\e-'^'\r') (3.6) 

(Af*^*^ corresponds to the one with Aif^"^ instead of Aii) and it satisfies the equation 

{ds-d^,+Vi{r)}Ai{r,r';s) = 0. (3.7) 

Since we are interested in the large / behavior of the Green function, let us rescale the 

Vi = l^Uu s = ^. (3.8) 

Then ( 13. 7p becomes 

i^dt-^d^+Ui{r)^Ai(r,r';^^=0. (3.9) 

From this equation one may readily recognize that the large / expansion is of the same nature 
as the derivative expansion which is an expansion in the number of derivatives on V/. 



In the rest of this section, we develop the derivative expansion of the above proper-time 
radial Green function. The resulting series will be identical with the large I series considered 
in [lO|, but a simpler derivation is given here. Introducing the momentum variable p, the 
proper-time Green function can be cast into the form 

A(r,r';s) = H ^{r\e-'^'\p){p\r') 

J-oo ZTl 

= r ^g-^[-9.2+V,(r-)]g-ip(r-r')_ (^3 ^Q) 

J -00 277 

After moving the last Fourier factor e"*^^''"'''-' to the left of the differential operator dr, one 
may set the coincident limit r' = r and get 

00 111 
dp 



r P-sv-k{t,p-s), (3.11) 

J ~oo ZTT 



where a new function K{r,p;s) is introduced. K{r,p;s) satisfies the following differential 
equation 

[ds - dl + 2ipdr + V/(r)} i^(r,p; s) = (3.12) 

and the boundary condition K{r,p; 0) = 1. 

Now, we introduce an auxiliary expansion parameter A in 03.12p . i.e., consider 

(ds - \^dl + 2iXpdr + Vi{r)]K{r,p- s) = 0. (3.13) 

Then, taking A as an expansion parameter, a series solution to the above equation can be 
found with the ansatz: 

K{r,p; s) = e"'^'^'^) [l + Xbi{r,p; s) + X%ir,p; s) + ■ ■ ■] . (3.14) 

Plugging (I3.14P into fl3.13p . we find the recurrence relations 

dsbk = -2ip {dr - SVI} + {d^. - 2sVidr - SV/ + {Vif} 6fc-2, (k > 2) (3.15) 

together with = 1 and dsbi = 2ipsV'i{r). With the boundary conditions 6fe(r, p;0) = 0, 
these recurrence relations may be used to determine the coefficient functions 6^ (A; = 1, 2, ■ ■ ■). 
Some leading terms satisfying fl3.15p are easily found. Note that, bk{r,p;s) is a simple 
odd/even polynomial of p when k is odd/even and hence it vanishes as we integrate over 
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p for all odd numbers of k. Here we report a few leading b^s with even numbers of k 
(A; = 2,4,6): 



62 = 
h-- 
+ 
+ 



WY + 



2p's 



,2 o3 



V,' 



I 5 



(3.16) 



' 4 8 

p S 



p^s^ 



' 4 

p S 



24 6 
15 



,2.6 



3 

2/s5 
15 



+ 



+ 



p^s^ 

13pS^° 

1^^6 48~~ 
'2p^s^ 124pS^ 



5£^ 

12 

c3 



(4) 

2S7p^s^ 
540 



3 

'2pS6 
9 



36 

13^2^5 

15 



13s5^ 



/n2 



+ 



24 



(Vi 



.//\2 



+ 



(3.17) 



7^ 
60 



vr(v;)^ + 



2^,5 



315 



4 6 
p S ^ p S 



9 



24, 



V, 



(6) 



. 45 

p6 ^8 



19s^\ 



315 
32pS^ 



+ 



43ph' 3.n (5) 
~60~~20r' '~ 



, 720 



p4^11 ^ ^2^10 



.9 



72 



36 



,2^,6 



+ 



17p^s 



60 ; 

82/s7 
21 



63 
287pS9 
270 

VM'' + 

_ 301s6' 
360 , 



(3)^2 



+ 



7p^s^ ^ 25p2s7 



139s^ 



493p^s« 25s' 
+ — + 



' 6 9 

p S 



240 ' 56 
109pS* 



^ 15 180 



+ 



473p2s^ 43s' 



15 



28 
254pS^ 



420 



180 



47pS9 ^ 359p2s8 



45 

271s7 



+ 



720 ^ 
269^2^6 



315 180 
/2/s9 31pS* 



90 



360 



1260 



9 15 
M'\vir. (3.18) 



From the series solution of K{r,p] s) we perform the Gaussian integrations over p. Then 
the proper time radial Green function at the same points, A(r, r; s), is determined as 



-sV, 



A(r, r; s) 
+A6 



A/ivrs 

06 



1 + A^I 

lls^ 



12 



/\2 



288^ 360 



6 

j4 

4"0 



(v;)X + li(v;')^ + |^v;vf) 



60 



10368 
43s6 



/n6 



8640^ ' 10080^ ' ' ^ 252^ ' 



61s6 
15120 



yi/M;(3) ^^n;'A2T,(4) 23s^ 



2520 ' ' ' 1008 



5040 



19s5 
2520 



+ A 



(3.19) 



280 ' ' 840 

Clearly, this A-series is organized according to the total number of derivative on V;. Of 
course, we can set A = 1 now. We also remark that, ignoring the last sixth order term 
in (13.191) . the quantum mechanical WKB series used in the approximate evaluation of the 
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instanton determinant in [17| can be obtained from fl3.19p if we substitute the relevant 
expression for the potential V{r). 



IV. LARGE L EXPANSION OF THE HIGH ANGULAR MOMENTUM PART 



In this section we use the derivative expansion (13.191) to derive the large L expansion in 
(12. 9p . To this end it is necessary to identify the structure of the renormalization countert- 
erms, 6Q, first. For simplicity we use the dimensional renormalization method. Setting the 
dimension of space-time to be (i — 2e, a dimensionally regularized expression of the functional 
determinant is, using the proper-time representation, 



°° ds e'^-^^f/^ 



„ oo 

J rf"'x^afc(x,x)s^ 



(4.1) 



s (47r)'^/2s^/2- . 

where n, introduced for dimensional reason, carries the dimension of mass. In (14.11) . afc(x, x) 
are the heat-kernel coefficients. Even if many of them are known in explicit 



(fc = l,2,3. 



forms, two leading coefficients, ai 
The integration over s yields 



V and a2 = ^V"^ — ^d'^V , are sufficient for our purpose. 



m 



k=l 



m 



2k 



J d'^x afe(x,x). 



(4.2) 



When the dimension of space is an odd number, above expression is finite at e = and does 
not require any counterterm. However it has a pole when the associated dimension is even 
and, particularly for d = 2,4, it has the structure: 

1 rd^x , 



47r 
d^x 



ai(^x,x), id = 2), 
[-m^ai(x, x) + a2(x, x)], (rf = 4), 



(4.3) 
(4.4) 



e J IGvr^ 

in the e — >■ limit. For d = 2 and 4 we here choose the renormalization counterterms, 
assuming the minimal subtraction scheme, as follows: 

1 n 



6n = —(^--7^) y'rf2xai(x,x), id = 2), 

m^ai(x, x) + a2(x, x) , (rf = 4), 



sn 



1 



--lE 



where 'Je is Euler's constant. With these counterterms, fin has a finite expression 

{H^s)'e-"' ' I dr 

--L+1 



l^H = liin 



/ _ (^2^)^e-™ ' dr Y: 9i{d) {A,(r,r; s) - Af-(r,r; s)} + 5^1 
Jo s Jo ^ J 



(4.5) 
(4.6) 

.(4.7) 
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From the expression (14. 7p and the derivative expansion of the proper-time radial Green 
function Ai{r,r] s) obtained in the previous section, the desired large L expansion in fl2.9p 
can be derived. After plugging fl3.19p into fl4.7l) . we may perform the / summation first. 
Note that the sum has the structure 

oo 

J2 e-^(^2i2+^i'+^o)(co + Ci/ + C2/2 + ■■■). (4.8) 

l=L+l 

This kind of summation cannot be done explicitly; but, with the help of Euler-Maclaurin 
method, it is possible to obtain the large-L series expansion. In computing the asymptotic 
expansion of this sum, the most useful form of Euler-Maclaurin formula is 

°° /-oo 1 1 1 , . 

^Ey(0 = I fm-^f{L)--nL) + —J<'\L) + ..., (4.9) 

assuming that /(/) and its derivatives vanish at / oo. When /(/) is of the form fl4.8l) . the 
integral in Euler-Maclaurin formula can be performed in terms of the exponential function 
and the error function along with a polynomial of L. Regarding s to be of an order of 
we find that taking a derivative of f{L) increases the power of Therefore the rest of the 
series in Euler-Maclaurin formula is an asymptotic large-L series. See Appendix C in [l^ 
for more details. After using Euler-Maclaurin formula and changing the integration variable 
s with -p-, we can find the large L series for the high angular momentum contribution in the 
following form: 

oo oo 

Y: giid){Mr,r;s)-A'rir,r;s)}= P-nL-""- (4.10) 

l=L+l n=2-d 

Because of the presence of the degeneracy factor, explicit forms of P-nS and further evalu- 
ation depend on the dimension of the space. However, since the procedure itself is basically 
the same, we will present the related calculation in detail for d = 2 and only the final results 
for (i = 3, 4, 5 below. 



A. 2D 



Note that the degeneracy factor is simply gi{2) = 2, which is independent from /. The 
large L expansion in (14.101) starts from Pq. Some of the leading coefficient functions P_„ are 
explicitly evaluated as 



Po 



-erfc 



(4.11) 
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Other coefficient functions have similar structures: one part being e"*/''^ times a polynomial 
in V or its derivatives, and the other part erfc(^) times a different polynomial in V or its 
derivatives. One may perform t integration with these explicit forms of P-nS. In performing 
this process with the first term, Pq, we find a pole term: explicitly. 



/ 

Jo 



with 



« = Wl + -^. (4.15) 



-ee-'T^ = ^r(n), (4.16) 
W#erfc[^] ^ 2(0'"^i^,Fi(n,n + l/2;n + l;l-K^), (4.17) 



This i divergence is canceled by the renormalization counterterm 5Q. In the evaluation of 
other terms on the other hand, no such divergence arises and so the limit e — > can be 
taken safely. In the evaluation of these terms, following integral formulas are useful: 

Jo t \r I \2j r(n + 1) 

where 2P1 is the hypergeometric function. Here note that we do not expand the function 
u — yr+^p^ as a power series of (even for a large value of L), for this kind of 
expansion breaks down when m or/and r get large. Keeping the function u as a whole, we 
can maintain the uniform nature of our large L expansion. 

After the Z-sum and the t integration, we can generate a j series for the large partial 
wave contribution fl^ as desired. The leading term of this series is 

which is the only nonvanishing term as we let L ^ 00 (since Qo — 0). Other terms vanish 
in the limit L — > 00 like {n >1), but they can be important for L not too large. Some 
of those secondary leading terms, needed for the fast evaluation of Jlren, are found to have 
following forms: 

= ~rV, (4.19) 
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g_2 = r l-er^wV^ + 2rV\/" + Ir^iu^ + m + 1)m V 

24'U'^('U + 1) 

-r(M^ + M2-3M-3)r}, 
g„3 = ^ fer^M^V^ _ ^r'u^V" - Qr\^V' - ?,riu^ - Qu^ + 5)^} , 

48?;/ I- J 



(4.20) 
(4.21) 



Q- 



— |80r^(2u + l)uV^ - 40r^(2u + l)u^{V') 

1)^ L 



1920L%ii(u + 

+8r^(2M + I)mV(^) + 16r^(2M3 + 4^^ + 6m + 3)mV(=^) + 60r2(M + lf{u^ - 5)mV2 
-r(M + 1)^(81^6 - 1185n^ + 2695^^ - 1575)^ - 80r^(2M3 + Au^ + Qu + ?,)u^VV' 
+4r2(4M^ + 18m^ + 32m^ - 139m^ - 310m^ + 20^^ + 350u + 175)uV 
-4r^(4u^ + 13m^ + 22m^ - 44n^ - llOn - 55)u V" - 80r^(2M + l)u^VV"] . (4.22) 



The formula given in [l^ can be reproduced immediately, utilizing only the piece Q\og 
above. First note that their result in 2D can be written in the form 



f^DK = lim 

L^oo 



Y.2ni + / drrVl-Y,- 

.1=0 I 1=1 ' 



In 



fir 
Y 



lE 



(4.23) 



Now, using the relation 



-|:^.n(f)+..^..>(|^)^o(^-'). 



(4.24) 



one can see that the second part of fl4.23p indeed corresponds to Qiog- (Here observe that 
M ^ 1 as L ^ oo). This clearly shows that the result of Dunne and Kirsten, ^dk, is equal 
to 



Odk = liiii 

L— >oo 



^ /-oo 



(4.25) 



i.e., our expression without any -^--suppressed terms involving Q^i,Q_2, 



B. 3D 

In 3 dimensions no ultraviolet divergence problem arises in the dimensional regularization 
procedure and hence no renormalization counterterm is necessary, i.e., 5VL = 0. Now the 
degeneracy factor is (7/(3) = 2/ + 1 - it grows linearly with I. Hence the series in (14.101) starts 
from Pi. Using the integral formulas in (14.161) and (14.171) . we can perform the t integration 
explicitly to obtain following expressions for the Q's: 

= -urV, (4.26) 
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--rV, 
u 



(4.27) 



g„i = {-Gr^uV^ + 2r3uV" + 2r'^{2u'^ + l)uV - 3r(4u^ - Sm^ - 1)^} , (4.28) 
Q_2 = {erVV^ - 2rVl/" - 6rW + 3r(2M^ + Sm^ - h)v] , (4.29) 

+16r^(2u2 + 3)mV(^) - 20r3(4M^ - Qu^ - ll)u^V" - 20r2(30«^ - lOn^ - 35)uV 
-SQr^u^VV" - 80r^(2u2 + 3)mV\/' + 15r(88M'^ - 235n^ + 42^^ + 105)\/} , (4.30) 

g_4 = (l20rV(V')^ - 60r3(6M^ + ISm^ - ?,h)u^V^ - 24rV\/(^) - 240r5MV^ 

-240rVl/(^) + 20r^(6M^ + 39m^ - 77)mV" - mr'^{Au^ + 18m^ - W^u^ + 105)mV 
+240rVl/V" + 1200rVl/V' + 15r(48M^ - 20m'^ - 1015m^ + 2142m2 - 1155)^} . (4.31) 

These explicit results can be utilized for fast evaluations of i^ren in = 3. Note that there 
is no In L related term in this dimension. 

If one does not care much about the fast convergence of the expression (in the limit 
L — > oo), the renormalized quantity fircn can be found using only the terms Qi and Qq 
above, i.e.. 



lim 

L— >oo 



^ poo 

^2(21 + 1)^1+ / dr{Q,L + Qo) . 
i=o J-, 



On the other hand, the 3D formula for the same quantity given in iJ] reads 



(4.32) 



^DK = 

L— >oo 



^ roo 

V(2/ + l)fii-/ drrV{r){L + l] 
1=0 



Observing that u = 1 + 0{L ^), one may easily see that 

{QiL + Qo) = -rV{r){L + 1) + 0{L-': 



(4.33) 



(4.34) 



and therefore two quantities in fl4.33p and (14.321) are identical. 



C. 4D 

The degeneracy factor in 4 dimensions is gi{4^) = (/ + 1)^ - it grows like a quadratic power 
in I. The large L series in (14.101) starts from P2 in this case. The integral formulas in (I4.16P 
and fl4.17p enable us to perform the t integration again. After some calculations we have 
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found the following results for the unsuppressed quantities: 

gi„g = -1 In (-J^) (6mW + 3rV2 - SrV - r V"), (4.35) 
24 y (it + Ij-iv y 

Q2 = -liu^ + 2u-l)rV, (4.36) 

o 

= -IrV, (4.37) 
Qo = {-6rVl/2 + 2rVl/" + 2r2(3M2 + 1)mV' + r(-52M^ + 25u^ + 3)1^} . (4.38) 

Subleading terms in the j asymptotic expansion can also be identified, with the results 

g_i = {Gr^u^V'^ - 2rVy" - 6r^'uV - r{16u^ - 37u^ + Gu^ + 15)v\ , (4.39) 

LI 

g_. = ,,3,o^n(^. _ {240r^(K3 - l^V' - UOr^u' - l)u\V'f 

+2Ar^{u^ - I)mV(^) + UAr^{u^ - l)u^V^^^ + 60r^(52w^ - 127m^ + mu^ + 15)itV^ 

+4r3(18M^ - 260m^ + 635m^ - 22Su^ - 165)mV" - 720r^(M^ - l)u^VV' 

-l2r'^{<ou^ + blQu^ - lOllit^ + 320^^ + 175)wV - 2AQr^{u^ - l)u^VV" 

+3r(ii^ - l)^(3288w^ + 1605m^ - 7630i/^ - lb7b)v] , (4.40) 

g_3 = — ^ |-720r^MV2 + 360r^M«(V^')' - 72r^«V(^) - 720r^MV(2) 
7680ii^'* 

+60r3(16ii^ - lllit^ + 30it^ + 105)itV^ - 20r=^(16it^ - lllit^ - 421^^ + 231)iiV' 
-60r2(90w^ - 201w^ - 182^2 + 315)mV + 12Qr^u^VV" + 3600rVyy' 
+ 15r(96M^° + 63m^ - 2980m^ + 3010m^ + 3276^^ - 3465)1^} , (4.41) 
Q_4 = _ {r\2u^ - + [-1680^^ + 3360y(y')^ - 1008(y")^ 

+3360l^V" - ISAAV'V^^^ - &12VV^*^ + A^V^^^] + r^{2u^ - + 5)m^° [432^^^^ 
-7392yV" - 4032y\/(^) + 10080^ V] - 1680r^(w2 - lf{b2u^ - I2bu^ - 35)itV^ 
+5r(w^ - 1)^(64692m^° - 748223^^ - 1201788^^ + 7638246m'' - 4708704?/^ 

-1576575)V - 168r^(28'u^ - 260m^ + 1145m^ - 1083m^ - 355^^ + 525)m*(1^')^ 
-336r^(12it'^ - 2mu^ + 1145ii^ - 1167ii^ - llbu^ + 3^b)u^VV" 
+24r^(36ii'^ - 364ii^ + 1603i/^ - 1365m^ - 785^^ + 875)iiV(^) 
+ 1008r^(4M^^ + 850m^ - 2869m^ + 2645m^ - l^W - 525)mVV^' 
-144r^(8it^^ + 1190ii^ - 3843ii^ + 3030ii^ + mhv^ - 1050)itV(^) 
-A2r^{u^ - l)2(9864ii^ - 8415ii^ - 64645it^ + 544951^^ + 17325)itV^ 
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+2r^(432M^^ + 69048m^^ - 25641m^° - 1242409u^ + 2941010u^ - 1987510m^ 
-167265m2 + 412335)mV" - 6r^{UAu^^ + 46032m^^ - 228564m^^ - 659897m^° 
+4287475m^ - 6809390u^ + 3736110u^ + 153615^2 - 525525)mV} . (4.42) 

These can be used for fast evaluations of Qj-en in d = 4. For the comparison with the 
result of y^], we here give the formula presented in [l^ : 



Qbk = lim 

L— >oo 



.1=0 I 



/o°° drrV{r) /o°° drr^V{V + 2m^) 



2(/+l) 8(/ + l)3 

drr^ViV + 2m2) jln (^^^ +7^ + 1 

The I sum of the second term in the above equation yields 

rV{r) 1 



(4.43) 



(4.44) 



while the / sum of the next term, if combined with the last term, produces 



V{V + 2m^ 



.1=0 



fir 



V{V + 2m^ 



(4.45) 



On the other hand, our formula for Qj-^n using only unsuppressed terms (given in (HI 
Km ) is 



lim 



11=0 



^(/ + Ifni + / c/r(Qiog + Q2L'' + QiL + Qo) 



(4.46) 



+ Qs^^ + QiL + Qo 



\2LJ 1 8 24 I ^ M 4 4 2 



8 24^ ' ^ 



(4.47) 



The difference between Q^^n and (14.431) is 

f2ren ~ ^DK = 



"rfr-IrVlnf^UrV^ = 0, 



24 



2L 



(4.48) 



showing that our result is consistent with the 4D formula given in 14 1. 
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D. 5D 



The degeneracy factor in the 5 dimensional space is gi{5) = |(/ + + 2) (2/ + 3). The 
first term in the large L series in fl4.10p will be the term P3 now. The integral formulas 
in fl4.16p and fl4.17p are used in the t integration again. After the I summation and the t 
integrations, we have found the following results: 

= ^u{2u' - 3)rV, (4.49) 
Q2 = -^rV, (4.50) 



Qi = {6r%2u'^ - l)u^V^ + 2r^{l - 2u^)u^V" + 2r'^{-Su^ + Au^ + 1)mV 

+r(-138M^ + 47^2 + 3)1^} , (4.51) 

Qo = (er^M^V^ _ 2r3MV" - Gr^wV - 3r(28M'' - 34m^ + 9^^ + h)v] , (4.52) 
72m' J 

Q_i = \80r'(2u^ + lyv^ - A0r%2u^ + iy(V'? + 8r%2u^ + I)mV(^) 

11520u^^ 

+16r^(8M^ + 4^2 + 3)mV(3) + Q0r^{A6u^ - 47^^ - 5)u^V'^ 

-5r(1152M^° - 5600m^ + 5010m^ + 2469m^ - 2716m^ - 315)^ 

-80r^(8M^ + Au^ + 3)u^VV' + 4r^(32M^ - 214m^ + 247m^ + 55)mV" 

-4r2(32M* + 16u^ + 1172m^ - 1045^^ - 175)^2^ - 80r%2u^ + 1)^^^"} , (4.53) 

Q_2 = ^ ^ |-240r^MV^ + 120rV(l^')^ - 24r^uV(4) - 240rVy(^) 
5760m^'^ 

+60r^(28M^ - 102m^ + 45^^ + 35)mV^ - 20r^(28M6 - 102m^ + 21^^ + 77)u^V" 

-180r2(32M^ - 74m^ + 7^^ + 35)mV + 2A0r^u^VV" + 1200rSVV"' 

+ 15r(704Mi° - 1656m* - 1260m^ + 3745m^ - 378^^ - 1155)v} . (4.54) 

One may use these results for fast evaluation of firen in c? = 5. 



V. THE PREFACTOR IN FALSE VACUUM DECAY RATE 

We will here illustrate how the formulas obtained in previous sections can be used to 
improve the rate of convergence in the calculation of the prefactor in the false vacuum 
decay. We consider a simple four- dimensional scalar field theory described by the Euclidean 
action 

5ciM = /rf'x 1^(9^0)2 + f/(0)} , (5.1) 
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with 

f/(0) = ^(02-a^)'-^(0-a) . (5.2) 

Here the parameter e represents a constant external source, which serves to break the degen- 
eracy of the double-well-type potential. The potential U (0) has two nondegenerate classical 
minima, 0_ and 0+(> </>-), with U{(j)J) > U{(f)^). After expanding the field (f) about the 
false vacuum 0_ 

= 0- + y^, (5.3) 
it is convenient to rescale the field (p and the spacetime coordinates as 

2 

Tfl 

X = mx ; ip = — $ (5.4) 

in the dimensionless form. Here the parameters m and i] are related to the original couplings 
by 

m' = ^(3<pl-a') ; V = ^\<P-\- (5.5) 
Then the classical action in terms of these dimensionless quantities is 



(5.6) 



with the dimensionless quartic coupling constant a = 

The bounce ^cii'r), which determines the decay of false vacuum, is a solution to the 
nonlinear ordinary differential equation 

- <^'ci - ^<fd + <fci - I'^l + f $^1 = 0, (5.7) 

satisfying the boundary conditions 

$;i(0) = 0, $ei(oo) = 0. (5.8) 

It is hard to solve this equation analytically but one can always find a numerical solution. 
The false vacuum decay rate is denoted by 7decav and in the one-loop approximation it is 

n 

given by [I!?] 

7decay-|^ 27r ) \det[-d^ + U" ^)] ) ' ^ ' 

where the prime on the determinant means that the zero modes (corresponding to transla- 
tional moves) are removed. In the exponent the first term 5'[$ci] denotes the classical action 
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of the bounce and SQ denotes the renormahzation counterterms. Note that this quantity 
7decay iuvolves the functional determinant and thus it can be evaluated using the methods 
developed in this paper. The potential V{r) in (12.1 p is now fixed as 

V{r) = f/"($d) = -3<l>ei(r) + ^<l>Ur), (5.10) 

and so the effective potential Vi in the radial operator Aii for a partial wave with / becomes 

= ^^^^^1^^^^ - 3$ei(r) + ^$i(r). (5.11) 



For the comparison with the result of [1^, we will call the logarithm of the prefactor in 
7decay (with the oppositc sign) the effective action F. Then the partial wave expression for 
the renormalized effective action is 



^ ^ 1=2 



+ ldr I^Qiog + ilQ-n^^ + O (-1) (5.12) 

with Qi = lndet(7W; + m^)/lndet(7Wf'''' + m^). In this expression, Qq has a negative sign 
and its absolute value is taken. The contribution from the sector with I = 1 involves the 
zero modes related to translational invariance and it is, having been removed from the sum, 
written down separately. The analytic expression for that contribution has been found in 
[l^ . The factor | in front of (15.121) is introduced since we are considering a real single scalar 
field and (/ + 1)^ denotes the degeneracy factor. In the last term, $o = *&ci(0) and $oo is the 
coefficient of Ki{r)/r {Ki{r) denotes the modified Bessel function) in the large r behavior 
of ^ci{i^)- The coefficient functions Q's are explicitly given in our subsection IIV Ci 

In the case a = 0.5 we plot in FIG.[T]the numerical values for the right hand side of (I5.12p 
as a function of L = 1, 2, ■ ■ ■ , 100, taking -^-suppressed parts of the asymptotic expansion 
successively. The lowest plots in this figure is the case when we ignore all suppressed terms 
in the large-L expansion of the high angular momentum part. It shows the existence of the 
L — >■ cx) limit. But, since the rate of convergence is quite slow, it is difficult to find the actual 
limit value which is the desired value for the effective action F. Other plots represent the 
values obtained after incorporating -p-, and corrections successively and they show 
remarkable improvements in convergence as we include these corrections. Magnified form of 
this figure is shown in the FIG [21 A relatively small number of L, for instance L ~ 30 in 
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FIG. 1: Plot for the sum of the low angular momentum (l < L) and the high momentum part 
when we take L = 1, 2, • • • , 100. The (brown) pluses denote the values as we ignore all terms of 
0{j^) and beyond. Slow convergence is evident. Solving 100 differential equations to determine the 
determinant for each partial wave (for the case with L = 100) is not enough to approach the limit. 
The blue crosses denote those after incorporating 0{j^) corrections and show the marked improve- 
ment of convergence already. The (red) diamonds, (gray) squares, and (black) dots represent the 
cases obtained after we incorporate -p-, and corrections successively. 



the a = 0.5 case, produces a good convergence and thus gives us a very accurate number of 
the effective action when / (ir[X]^=i Q-nL'""] is added. 
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FIG. 2: The same plots with Fig [T] but the scale is magnified by 250. We can see only the (red) 
diamonds, (gray)squares, and (black) dots which represent the values after incorporating -p-, -p-, 
and corrections successively. One may conclude that L ~ 30 is enough to get the limiting value 
with a very good precision. 

VI. RADIAL OPERATORS IN GAUGE THEORIES 



In this section we consider an application of the large L expansion in (12.91) to gauge 
theory. In p^] it was used to calculate the renormalized effective actions in classes of 
radially symmetric background gauge fields. We present here the explicit forms of some 
coefficients functions Q-n which were announced in that paper. 

The one-loop effective action in a gauge theory is expressed by the logarithm of the 
functional determinant, 



= In 



(6.1) 



det[— 9^ + m^] J ' 

where = D^D^ and = d^ — iA^ is the covariant derivative operator. We consider two 
cases with SU(2) backgrounds gauge fields of the form 



(Case 1): A^(x) = 2'q^^aXuf{,r)—, 
(Case 2): A^(x) = 2{7]^^iUi)x^g{r)—, 



(6.2) 
(6.3) 
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where iJi^v = 1,2,3,4, rj^j^ya (or rjfj_yi) are the 't Hooft symbols and Ui a unit 3-vector. 
The functions, /(r) and g{r) are unspecified so that general (radially symmetric) back- 
ground fields can be studied. Following we define the angular momentum operators 
La = —-^Vfj-vaX^di, and the SU(2) isospin generators = which satisfy the commuta- 
tion relations [La,Lb] = ieabcLc and [Ta,Tb] = ieabcTc- These operators carry the quantum 
numbers characterized by / (= 0, |, 1, ■ ■ ■) and t (= |). Note that the angular momentum 
operators defined here are different from those of previous sections so that the quantum 
number / can take half-integer values as well as integer values. (The quantum number / in 
previous sections corresponds to 21 in this section). Because of the radial symmetry, as in 
previous sections, the one-loop effective action can be written as a sum over one-dimensional 
radial determinants. 



A. Case 1 

The form of the gauge field (16. 2p is inspired by the instanton solution (which corresponds 
to /(r) = ^2^p2 with the size parameter p) and it carries genuine non-Abelian nature. The 
covariant Laplacian —D^ involves the isospin-orbit coupling term and then {J^ = La+Ta), 

and = | are conserved quantities. Therefore, partial waves are specified by the 
quantum numbers (/, j). For each sector, the radial differential operator associated 
with the covariant Laplacian —D^ assumes the form 

Mii,) = -d'r+Vii,,) (6.4) 

with the effective potential 

V(,,) = ^ ^ + 4/(r) + 1) - /(/ + 1) - ^} + 3rV(r)^ (6.5) 
The corresponding radial operator for the free Laplacian —d"^ is 

M'r = -dl + Vr = -dl + ^^^^^^f^^^^ (6.6) 
Introducing the one-dimensional radial determinant ^{ij) 



^ J^^B^^ll^]^ (6.7) 



22 



the low angular momentum part of the one-loop effective action can be written as 

Ol= + + + (6.8) 

and the corresponding high angular momentum part as 

oo 

f^H = E (2^ + 1)(2^ + 2) + + (6.9) 

with the renormalization counterterm 



sn 



--lE 

e 



dr — 



4rV(^)^ - 8r^f{rf + Sf{rf + 4r/'(r)/(r) + r'^firf] . (6.10) 



[The Pauli- Villas regularization was employed in 15| but we have changed it to the dimen- 
sional regularization scheme in this work]. In fl6.8p and (16.91) we have rearranged j(= / ± |)- 
sum by combining two sectors of (/, j = I + \) and (/ + \,i = I) with the same degeneracy 
factor {21 + 1)(2/ + 2). See (lol for details. 

Using the WKB series described in Sec JIIIt the high partial-wave contribution (16. 9p can 
be calculated analytically in the form of large-L asymptotic series. The calculational step is 
almost the same as that described in Sec. IIVI A part of the result was presented in |l^ by 
the form 



/•oo f °° 1 

^^H = / dr\ Qiog In L + ^ Q-„i.-" • 

[ n=-2 J 



(6.11) 



The first few terms, i.e., Qiog, Q2, Qi 



and they are already presented in 
1 



, Qo and Q-i are given by relatively short expressions 



15|. We here report the explicit expressions for Q_2 and 



Q- 



— {-80(16m^ + 27u^ - 81m^ + A5u'^ - 7)m^G^ 



7680^13(^2 

-f20(32M^^ + 792m^° - 5175m* + 8910m^ - 5140m^ + 266^^ + 315)u^G^ 
+15(m^ - 1)2(3472m® - 5619m^ - 7399m^ + 7623m^ + 1155)6" 

-40r(12M9 + 54u^ - 171m^ + 140^^ - 35)^^(6^)' - 60r^{4u'^ + 6u^ - 15u'^ + 5)u%G'Y 
+40r2(4u^ + 9u^ - 18^^ + 5)u\G^y + A8r^u%u^ - l){Au^ - 5)0^'^^ 
+4r(1760M^° - 12494m* + 25779m^ - 19595m^ + 2975^^ + 1575)^^6" 
-4r2(440M* - 1889m^ + 2019u^ - 185^^ - 385)n^G"' - 24:r\u^ - l)u^G^^^ 
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+2Ar\u^ - l)u^{2>{H"f + AH'H^^^) - IGr'^iu^ - 36u^ + 35)u'^H'H" 
+Ar^{6u^ - 930m^ + 1854m^ - 685m^ - 2A5)u\H'f} , (6.12) 
Q_3 = ^ ^ 1-720(9^6 - 45u^ + 35^^ - 7)u'^G^ - 60r\l8u'^ - 75u^ + 35)u%G'? 
+60(72m1° - 1431m*^ + 3170m^ - 420m^ - 2562^^ + 1155)u^G^ 
+15(528m^2 - 3833m^° - 8245m® + 41454m^ - 29442m^ - 15477m^ + 15015)G 
-40ru\lQ2u^ - 855m^ + 980^^ - 315) (G^)' + 40rV(27M^ - 90^^ + 35) (G^)" 
+12r(160Mi° - 3534m® + 10165u^ - 4655u^ - 7875^^ + 5775)^^6" 
-12r2(40M® - 449m^ + 165m^ + 1435m^ - 1155)u^G"' + 48r^u^(12u^ - 45m^ + 35)G'(^) 
-24r^(3M2 - 5)m®G(^) - 120r^M®(3(if")' + AE'R^"^^) + SOr'^u^SGu^ - A9)H'H" 
-20r2(174M^ - 414m^ - IGlw^ + AAl)u^{H')'^} , (6.13) 



where u{r) = y 1 + H{r) = r^f{r) and G(r) = H{r){H{r) - 1). The Q_4 term was 
also used for the evaluation in |15], but the expression for Q_4 is quite long and its actual 
numerical value is rather small in most cases. So we do not present it here. 

In 15| it was shown that incorporating the combination Q-n-^ " made the summa- 



tion over / converge dramatically fast. For instance, FIG. 4b in 15| clearly demonstrated 



the changes when each of Q-i, • • • , <5-4 terms was added one by one. 



B. Case 2 

The second case is quasi-Abelian. The field has a fixed color direction and only 
{LF'^L^^T^} are conserved quantities. Partial waves for the covariant Laplacian —D^ are 
classified with the quantum numbers {l,l^,t^) (for = ^3 = ±|). The radial 

operator for a given partial wave has the form: 

Mm,,t,) = -a^ + V(i,z3,t3), (6.14) 

V(M3,t3) = ^ ^ + 89{r%h + r'9{rf. (6.15) 
The one-loop effective action in this case can be written as 

L I cx) I 

^= E E E (2^ + WM3,t3) + E E E (2/+ WM3,t3), (6.16) 

Z=0,i,l,- ^3 = -' t3=±i l=L+^ ^3 = -' t3=±i 
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where ^In^^^t^) is the one-dimensional radial determinant defined by 



/ detl- Df,, , ^ + m'^] , 



As in Case 1, the result of renormalized large partial wave contribution can be written 



in the form (16. lip and (^2, ■ ■ ■ , Q-i as well as Qiog terms were given in [15|. The expressions 

for Q_2 and Qs are 

= 2304oj'(^u^ - 1) {-240«'(«' - ^rrV + 240u%u' - 1)(4«2 - 3)r^/(^')' 
+mu''{u^ - l)V^(7V + 480m^(m' - l){Qu^ + 3u^ - l^g^g' - ASgu^'iu^ - l)/(?(^) 
+60m2(52m^° - 303m^ + 678u^ - 560m^ + 70u^ + 63)r^^^ 

+24m^(m^ - 2m^ + l)r^(3(/)2 + Agg^^^) + 16u%7W - 72u^ - 54^^ + 55)g'g"r^ 
+4m^(442m^ - 520m^ + 1606m*^ - 2043m^ - IOm^ + 525)r^(^')^ 
+8m^(216m^ - 260u^ + 1043m^ - 1599m^ + 215^^ + 385)r^^/ 
+40m2(72m^1 - 208m^° + 202m^ + 1131m*^ - 1855m^ + 343^^ + 315)r^^' 
+5(2248m^2 - 12421m^° + 11733m^ + 20070m^ - 32214m^ + TllQw^ + 3A65)g^ 
+96m^(2m^ - 4m^ - + h)r^gg^^^ } , (6. 18) 

= 3072L15 - 3)(«' - ^f^'s' + 2AQu\l2u' - 35u' + 21^ g^g'f 

+480m''(3m^ - IOm^ + 7)r^g^g" + 480m^(18m^ - Ihu^ - TOm^ + Q3)r^g^g' 
-48m^(3m2 - hy^gg^^'^ + 60m2(16m^° - 207m^ + 960m*^ - 1022m^ - 504^^ + 693)ry 
-2Au^{<ou^ - 5)r^(3(/)2 + Ag'g^^^) - 16m^(216m^ + 270^2 - 385)rV/ 
-96u''(12m^ + ISu^ - 35)r3^^(^) - 8m^(80m^ - 699u^ + 2595m^ + 1645^^ - MQ^y^gg" 
-4m^(160m^ - 1038m^ + 3465u^ + 3220^^ - 4725)r2(^')^ 
-40m2(64m^° - 174m^ - 1605m^ + 2905m^ + 2583^^ - 3465)r^^' 

+5(160m^2 - 4107m^° - 735m^ + 49770m^ - 44982m^ - 46431^2 + 45045)^^} . (6.19) 

Since the radial operator explicitly depends on Z3, we have to calculate the radial determinant 
for each /s: this makes the amount of calculation grow very fast, i.e., by quadratic powers 
of L as L becomes large. Thus the effect of acceleration, rendered by incorporating the 
J- suppressed terms, is greater than other cases. 
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VII. CONCLUSION 



We have here presented an efficient and precise method for the calculation of functional 
determinants with radially symmetric differential operators. This method involves the par- 
tial wave cutoff technique in which the infinite partial wave summation is separated into two 
sectors as in (12.61) . The first sector is evaluated with the radial Gel'fand-Yaglom method for 
each partial wave. We have developed the large-L asymptotic series for the second sector, 
i.e., for the high angular momentum part. Combining the first sector with the unsuppressed 
terms of the latter series, the renormalized sum can be found in the L ^ oo limit. Including 
the subsequent (i.e., -^-suppressed) terms in the series also, we can get a precise value for 
the functional determinant with the choice of relatively small L (which means less computa- 
tional work for low angular momentum part). Certainly this greatly improves the efficiency 
of calculation. That is, we can get an result with estimated errors of ~ if we ignore all 
the terms suppressed by in the large L asymptotic expansion (12. 9p . However, keeping 
the summation up to Q-N-pr, the estimated error rate will be reduced to ~ j^w+t- So with 
a suitable choice of L and A^, we can efficiently calculate the functional determinant to a 
desired accuracy with a relatively small number of L. 

A generalization of this work to fermi fields should be important. It can be done by 
converting the Dirac operator into the squared second order form, and then by applying the 
method developed in the present work. However it should be possible to develop a more 
direct fermionic approach, studying the coupled ffist order equations, along the line of the 
partial wave cutoff method. Developing a similar partial wave method, to evaluate quantum 
corrections to the masses of solitons like vortices and magnetic monopoles in gauge theories 
is certainly an interesting problem. 
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